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Abstract 

Using the reduced density matrix formalism the exciton scattering in light-harvesting 
systems of purple bacteria is calculated. The static disorder (fluctuations of the site 
energies) as well as the dynamic disorder (dissipation) is taken into account in this 
work. Circular aggregates with 18 pigments are studied to model the B850 ring of 
bacteriochlorophylls with LH2 complexes. It can be shown that the influence of dis- 
sipation may not be neglected in the simulation of the time-dependent anisotropy 
of fluorescence. Also an elliptical deformation of the ring could be essential. 
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1 Motivation 



Highly efficient light collection and excitation transfer to the reaction center 
initiates the energy conversion in photosynthesis. This process takes place in 
the so-called light-harvesting antenna network. Particularly, the ring subunit 
of the peripheral light-harvesting antenna (LH2) [1] of purple bacteria has been 
extensively studied [2]. The very symmetric arrangement with short distances 
between the pigments gave new impulses to the discussion about coherence 
in the exciton transfer and the exciton derealization in LH2. Both could be 
reduced by dynamic and static disorders. 

Time-dependent experiments made it possible to study the femtosecond dy- 
namics of the energy transfer and relaxation [3] . Kumble and Hochstrasser [4] 
have presented a time-domain analysis of the effects of the static disorder upon 
the dynamics of optical excitations. The interpretation of time-dependent ex- 
periments on the femtosecond time scale requires a theory which incorporates 
static and dynamic disorder. The aim of the present paper is to extend the in- 
vestigation by Kumble and Hochstrasser taking into account the simultaneous 
influence of static and dynamic disorders after impulse excitation. In addition 
we calculate the time-dependent anisotropy of fluorescence not only for the 
symmetric but also for elliptically distorted rings. 



2 Model 

The Hamiltonian 

H = H® x + H ph + H cx _ ph + H s + H c i (1) 

describes the transfer of a single exciton with a transfer integral J along the 
ideal ring (H® x ), the independent heat baths for each chromophore (H ph ), the 



2 



site-diagonal and linear interaction between the exciton and the bath (if ex _ p h), 
the static disorder (H s ) with Gaussian distribution (standard deviation A) and 
elliptical distortion of the ideal ring (H d ). 

Diagonalization of the Hamiltonian H® x of the ideal ring leads to eigenstates \k) 
and eigenenergies Ej, = —2 J cos k. For a symmetrical coplanar arrangement 
of site transition moments /2„ dipole-allowed transitions populate only the 
degenerate k = ±1 levels of the ideal ring. If the ring is elliptically distorted the 
eigenenergies of the Hamiltonian iJ° x + H c \ are no longer degenerate as is the 
case for the ideal ring. With static disorder of the site energies being present 
(A 7^ 0), the stationary states, i.e. the eigenstates \a) of the Hamiltonian 
i7g X + H s , correspond to mixtures of \k) and an excitation will prepare a 
superposition of the \k) states. 

The dipole strength jl a of state \a) of the ring with static disorder and the 
dipole strength jl a of state \a) of the ideal one read 

N N 

A*a = c nf i n, f^a = ^ C n^n- (2) 

n=l n=l 

The coefficients and c l n are the expansion coefficients of the eigenstates of 
the ideal and disordered rings in site representation. 

Kumble and Hochstrasser [4] concluded, that in the case of pump pulse exci- 
tation the dipole strength is simply redistributed among the exciton levels due 
to disorder. So the amplitudes of site excitations and the phase relationships 
in the initial state are necessarily identical to that of an equal superposition 
of k — ±1 excitons of the ideal ring. Thus, generally, the excitation with a 
pump pulse of sufficiently wide spectral bandwidth will always prepare the 
same initial state, irrespective of the actual eigenstates of the real ring. The 
nature of this initial state is entirely determined by the selection rules of the 
ring without static disorder. The initial condition for the density matrix by 
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pulse excitation with the polarization e x is given by (Eq. (la) in [5]): 

p a p{t = 0; ei) = ^(e x ■ jl a ){e x ■ ftp), (3) 



where A = T, a (e x ■ jl a ){e x ■ /*«). 



3 Anisotropy of fluorescence 

Kumble and Hochstrasser [4] calculated the usual time-dependent anisotropy 
of fluorescence 

r ( f \ _ (^xxjt)) - (Sxyit)) / j \ 

U (S xx (t)) + 2(S xy (t)) {) 

where, e.g., 

(S xy (t)) = (| E % ■ M(e y ■ A)CcLe-^f ). (5) 

a,l,n 

The indices a and / label the eigenstates of the virtual and disordered ring, 
respectively. The brackets () denote the ensemble average and the average over 
the direction of the laser pulses with fixed relative directions e x and e y . 

To include the dynamic disorder which contributes to dephasing of the initial 
wave packet and promotes thermalization of the dephased populations one has 
to work within the exciton density matrix formalism [6] instead of using only 
the exciton wave functions 

S xy {t) = J P xy {u,t)dw (6) 



where 



Pxy(w,t) = A^^p lv {t){e y ■ jli>)(e y ■ fii)[6(u - uj vo ) + 5(uj - u l0 )]. (7) 
i i' 
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4 Density matrix formalism for exciton transfer and relaxation 

Provided that the exciton dynamics is not very fast, its coupling to the bath 
weak, and except for the initial time-interval t ~ ta {t d = dephasing time of 
the bath), the adequate equation for a factorized initial state is the Redfield 
equation [7,8]: 

ij t p(t)= l -[H,p(t)]+1lp(t), (8) 

The Redfield relaxation superoperator 1Z describes the influence of the thermal 
bath on the dynamics of the exciton. 

Capek applied several different ways of obtaining of convolutional and convo- 
lutionless dynamical equations for the exciton density matrix in the site basis 
[9]. After Markovian approximation they have the following form 

l l p (t) = ^[H,p(t)]-6n P (t). (9) 

We proved [6] the equivalence of the Redfield theory without secular approx- 
imation with Capek's equations after Markovian approximation. 

5 Results 

In Kumble and Hochstrasser's modeling, the anisotropy of fluorescence of the 
ring LH2 subunit decreases from 0.7 to 0.3 — 0.35 and subsequently establishes 
a final value of 0.4. Kumble and Hochstrasser concluded that one needs static 
disorder of strength A ss 0.4 — 0.8J to reach a time decay below 100 fs . 

Our results for the time dependence (t = tK/J) of the anisotropy of fluores- 
cence Eq. (4) in the symmetrical ring are shown on Fig. 1 with static dis- 
order A/ J = 0.4. In addition we consider two strengths of dynamic disorder 
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jo = 0.2, 0.4 entering the spectral density J(u) = Q(uj)j -^je~ ul/u}c with cut-off 
frequency u c = 0.2 J [10]. Inclusion of dynamic disorder leads to faster decay 
of the anisotropy of fluorescence during the initial stage. Smaller strength of 
the static disorder A/ J than predicted by Kumble and Hochstrasser would be 
necessary to guarantee the decrease of the anisotropy of fluorescence during 
the first 100 fs. 

It was concluded in Ref. [4] based on measurements by Chachisvilis et al. [11] 
that the time decay of the anisotropy of fluorescence during the first dozens of 
fs is temperature independent in the case of LH2 subunits. Our calculation for 
the symmetrical ring show that such result can be obtained only for A > 0.8. 
But because the time resolution of the experiments in [11] was not too high 
this very restrictive statement about the strength of the static disorder can 
only be made with caution. We expect that some temperature dependence can 
be seen using experiments with shorter laser pulses. 

Recent results obtained by single molecule spectroscopy [12] can only be in- 
terpreted [13] admitting the presence of a C 2 distortion of the LH2 ring. It 
has, up to now, not been concluded whether such C 2 distortion of the LH2 
ring is present also in samples in vivo. 

We have made calculations which take into account the possible C 2 distortion 
of the ring using the model C by Matsushita et al. [13] of elliptical distortion 
of the ring in which the transfer integral J is cosine modulated. Local tran- 
sition dipole moments lie tangentially to the ellipse. Our results for the time 
dependence of the anisotropy of fluorescence (4) in the elliptically distorted 
ring are shown on Fig. 2 with the static disorder A/ J = 0.4 and the strength 
of the dynamic disorder j = 0.2 for three values of the elliptical deformation 
V 2 /J = 0.05, 0.1, 0.2. 

It is seen, that the inclusion of an elliptical deformation leads to faster decay of 
the anisotropy of fluorescence and diminishes the influence of the dynamic dis- 
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order. Our calculation for the elliptically distorted ring (V 2 /J = 0.2 [13]) show 
that the temperature independent time decay of the anisotropy of fluorescence 
during the first dozens of fs can be obtained even for A 0.4 . 
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Fig. 1. Time dependence of the anisotropy of fluorescence r(r) for the symmetrical 
ring with the static disorder A/J = 0.4. The influence of the dynamic disorder is 
displayed by curves for low (ii), (iv) and room temperature (iii), (v) for jo = 0.2 
(ii), (iii), and jo = 0.4 (iv), (v) compared to jo = 0.0 (i). 
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Fig. 2. Time dependence of the anisotropy of fluorescence r(r) for the elliptically 
distorted ring with the static disorder A/J = 0.4 and dynamic disorder jo = 0.2. 
The influence of the elliptical deformation displayed by dotted lines for V2/J = 0.05, 
dashed lines for V2/J = 0.1 and long dashed lines for V2/J = 0.2 compared to 
V2/J = 0.0 (solid lines) for low and room temperature. 
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